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ABSTRACT 



We report the negative results from a search for 6.7 GHz methanol masers 
in the nearby spiral galaxy M33. We observed 14 GMCs in the central 4 kpc 
of the Galaxy, and found 3a upper limits to the flux density of ~9 mJy in 
spectral channels having a velocity width of 0.069 kms -1 . By velocity shifting 
and combining the spectra from the positions observed, we obtain an effective 3a 
upper limit on the average emission of ~lmJy in a 0.25 kms -1 channel. These 
limits lie significantly below what we would expect based on our estimates of 
the methanol maser luminosity function in the Milky Way. The most likely 
explanation for the absence of detectable methanol masers appears to be the 
metallicity of M33, which is modestly less than that of the Milky Way. 

Subject headings: ISM: molecules - individual (methanol); ISM: masers; Galaxies 
(M33) 



1. Introduction 



6.7 GHz methanol masers are the second brightest maser transition ever observed, 
and are typically much brighter than OH masers. The properti e s of methanol mase rs in 
the Magellanic clouds ([Sinclair et al.l Il992l ; lEUingsen et al.l Il994l : iBeasley et al.l Il996l ) are 
consistent with those of our Galaxy, given appropriate consideration for different galactic 
properties such as metallicity. The SMC has an oxygen abundance 12 + log(0/H) = 7.96 
( iVermeij &: yan der Hulstl 120021). which is a factor ~ 6 smaller than that of GMCs in the 
Milky Way (IPeimbert. Storey, fc Torres-Peimbertlll993l ). 
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However, 6.7 GHz methanol masers have not been discovered in galaxies beyond the 
Magellanic clouds. In particular, surveys have shown that there is no analog to OH mega- 
masers at 6.7 GHz fellingsen et alill994j : IPhillips et aDll998l : barling et al1l2003h . It would 
be of interest to detect methanol masers in a Milky Way like galaxy. If the number of sources 
detected were large, one could derive the methanol maser luminosity function since all sources 
would be at the same distance. Further, one could look for correlations of methanol masers 
with giant molecular cloud masses, other types of masers, etc. The Arecibo radio telescope 
offers unequaled sensitivity for targeted surveys for methanol masers, and the nearest spiral 
galaxy which can be observed with this instrument is M33. 

It is difficult to develop an optimal strategy to search for extragalactic methanol masers 
since the luminosity function of methanol masers in our Galaxy is unknown. This is mostly 
due to difficulties in determining distances to methanol masers, which is compounded by the 
kinematic distance ambig uity. We have adopted the following approach. We selected sources 
from the general catalog (IPestalozzi. Minier. fc Booth! 120051 ). eliminating sources wihtin 10 
degrees of longitudes and 180 degrees since these suffer from large uncertainties in kinematic 
distance. Assuming that all masers are at their near kinematic distance, we then scaled the 
so urce flux density to t he distance of M33, which we take to b e 840 kpc, averaging the results 
of iLee et al.l (120021 ) and lFreedman. Wilson, fc Madord (Il99ll ). which gives a result consistent 
with the distance determined by Kim et al. ( 20021 ) . The resulting histogram of "expected" 
methanol maser flux densities is showin in Figure [T] 

We thereby obtained an estimate for the distribution of flux density of methanol masers 
in M33, which of course must be regarded with considerable skepticism, as the original Galac- 
tic sample has significant biases, not to mention the unknown differences in the properties 
of massive star formation in M33 and the Milky Way. It is also inherently pessimisitic in the 
sense that some of the Galactic masers are at their far kinematic distance, and hence are 
more luminous than we have assumed. With these caveats, we found a "high flux density 
tail" for our hypothetical methanol maser population in M33 that extends from 100 mJy 
down to 1 mJy. Approximately 16% of the total number of masers would be expected to 
have flux density greater than 1 mJy, with about 6% having flux density greater than 3 mJy. 
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Fig. 1. — Distribution of expected flux densities of methanol masers in M33 based on observed 
masers in the Milky Way. 
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Observations 



The search for methanol masers in M33 was carried out using the 305 m Arecibo tele- 
scope^and the C-band high receiver JPandian et alJbood ) between 2005 July 9 and July 21. 
The data were taken in standard Arecibo position switched mode in which the source is 
observed for five minutes, and then a reference position offset by 6 m in Right Ascension is 
observed for an equal length of time. The spectrometer used was the interim correlator. Two 
boards recorded the two orthogonal linear polarization data at 0.069 kms -1 resolution (3.125 
MHz bandwidth with 2048 lags) while the two other oboards recorded the same data at a 
higher velocity resolution of 0.034 kms _1 (1.5625 MH z bandwidth with 2048 lags). All the 
boards were centered on the CO velocity as given by lEngargiola et al.l (120031 ). The FWHM 
beam width of the telescope at 6.7 GHz is 40", which corresponds to a linear size of 160 
pc at a distance of 840 kpc. System temperatures varied between 23 K and 34 K, with ap- 
proximately 20 percent difference between the two polarizations. The data were converted 
to flux density using the elevation dependent gain curve for this frequency, with a typical 
conversion factor of 5 K/Jy. 

Based on the correlation between methanol masers, massive star formation, and Giant 
Molecular Clouds GMC) in the Mil ky Way, the posi t ions i n M33 observed were the most 
massive clouds given in Table 1 of lEngargiola et al.l (120031 ) . Each of the first 14 clouds, 
including all of the most massive GMCs in M33 (M > 4xl0 5 M ), was observed for a 
total of 30 minutes on plus 30 minutes for the reference position. The resulting rms for 
the low resolution data after averaging the two linear polarizations is 3 mJy, and for the 
high resolution data is a factor of y/2 higher. No spectral features having greater than 3a 
significance were observed. 



3. Data Analysis 



The observation of 14 massive GMCs in M33 has yielded no methanol maser detections 
at a 3a level of 9 mJy. Based on comparison with the expected distribution of flux densities 
based on 6.7 GHz methanol masers in the Milky Way, this is s urprising. There are estimated 
to be approximately 1200 methanol masers in the Milky Way (Ivan der Waltll2005l ). while the 
number of GMCs out to 8 kpc from the Galactic Center is 4000 to 5000 (Solomon & Sanders 
1980; Scoville & Sanders 1987). Having observed 14 GMCs in M33 with mass between 2 and 
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7 xlO 5 solar masses, if the relationship between methanol masers and GMCs were the same 
as in the Milky Way, we would have expected to have made several detections. 

Assuming that all the GMCs do exhibit methanol maser emission at a weak level, one 
can combine the data for all GMCs to get a better upper limit on the emission. This assumes 
that the emission occurs at the same velocity offset with respect to the velocity of the central 
channel in any given spectrum, which is not necessarily true. To allow for different velocity 
offsets for the location of emission in different GMCs, we follow the following procedure. 
We first select a section of the spectrum around the velocity of CO emission for each GMC 
(~ ±40 km s^ 1 ), which includes the largest offset one would expect between the methanol 
maser and the molecular line emission (expected to be less than 15 km s" 1 ). Each sub- 
spectrum is then cross-correlated with a Gaussian whose width represents the expected 
linewidth of the maser emission (~ 0.30 kms -1 ). Each sub-spectrum is then shifted by 
the offset between the central channel and the channel of peak cross-correlation. Following 
this procedure, the shifted data for all GMCs are added together. We denote the resulting 
spectrum as "A" , and is shown in Figure [2j The above procedure will create a line feature 
resembling the correlation template (a Gaussian signal in this case) from the constructive 
alignment of purely random noise. Hence, the procedure is repeated for a different portion 
of the spectrum, where one does not expect to see any methanol maser emission, and the 
resulting spectrum is denote "B" (Figure [3]). We then compare the A and B spectra. A weak 
emission feature would be manifest by the peak in the spectrum A being higher than the 
peak in spectrum B by a statistically significant amount. Since this is not the case as can be 
seen by comparing Figures [2] and [31 we conclude that there is no detectable methanol maser 
emission in the data. The resultant 3 a upper limit on the average emission is ~ 1 mJy 
This is a full order of magnitude below the flux density that we would expect to be common 
in M33 from the extrapolation of our admittedly uncertain results in the Milky Way. This 
difference, if real, would certainly suggest some large-scale difference in the molecular clouds 
or massive young stars in M33 compared to the Milky Way. A much more meaningful result 
could be obtained if we had better knowledge of the methanol maser luminosity function in 
our Galaxy. 



4. Discussion and Conclusions 

It is difficult to be highly quantiative about the lack of detection of methanol maser 
emission in M33 beyond the statements given above. Given that the 40" Arecibo beam 
subtends a region 160 pc in size at the distance of M33, we should be sensitive to a methanol 
maser located anywhere within the individual giant molecular clouds being observed. Based 
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Fig. 2. — The bottom panel shows the result of adding the spectra of 14 GMCs that have 
been shifted in velocity to the velocity of peak cross-correlation with a Gaussian template. 
This spectrum is denoted spectrum A. The x-axis of the plot shows the velocity offset from 
that of the CO emission line. The top panel shows the shifted individual GMC spectra, each 
of which is smoothed to a linewidth of 0.25 kms -1 . The side panel shows the amount of 
velocity shift applied for each spectrum. 
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Fig. 3. — The bottom panel shows the result of adding the spectra of 14 GMCs that have been 
shifted in velocity to the velocity of peak cross-correlation with a Gaussian template. This 
analysis has been done over a portion of the spectrum where no methanol maser emission 
is expected and the result is denoted spectrum B. The x-axis of the plot shows the velocity 
offset from that of the CO emission line. The top panel shows the shifted individual GMC 
spectra, each of which is smoothed to a linewidth of 0.25 kms -1 . The side panel shows the 
amount of velocity shift applied for each spectrum. 
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on our estimate of the Milky Way methanol maser luminosity function, we would expect a 
significant fraction of masers to have a flux density > 10 mJy, significantly higher than the 
individual limits we have obtained, and an order of magnitude greater than the averaged 
limit derived by combining the results from the 14 GMCs observed. 

Explanations of the rarity of methanol masers in external galaxies have focused on (1) 
insufficient methanol density over path where amplification cou ld take place and (2 ) insuf- 
ficient pumping to invert the methanol transition in question (IPhillips et al.lll998l ). Both 
of these can result from low metallicity. A reduced abundance of oxygen, for example, will 
likely reduce both the abuance of methanol (CH 3 OH) and of dust, which is required to con- 
vert the short wavelength radiation from massive young stars to the infrared wavelengths 
required for maser pumpi ng. The rarity of ma sers in the Magellanic clouds has been dis- 
cussed in similar terms by iBeasley et al.l (119961 ). The O/H ratio in the Magellanic clouds is 
dramatically lower than that of the Milky Way (see discussion in Beasley et al. 1996 and 
the more recent measurement by Vermeij & va n der Hulst 20 2). F or M33, the O/H ratio 
is only slightly less than that of the Milky Way. IVilchez et al.l (119881 ) determined that 12 + 
log(0/H) = 9.0 at the center of M33, falling t o ^8.5 at dista n ces be tween 2 and 5 kpc from 
the center of the galaxy. The best fit line of ICrockett et al.l (120061 ) to their new data plus 
previous data has a very small slope of -0.12 dex kpc -1 and 12 + log(0/H) = 8.3 at 5 kpc 
from the center of M33. These valu es may be compared to Galactic values of 1 2 + log(0/H) 
= 8.51 for Orion and 8.78 for M17 ( Peimbert. Storey, fc Torres-Peimbertlll993l ). The GMCs 
we have observed are located between 0.5 and 4 kpc from the center of M33, with a mean 
distance of 1.9 kpc. We conclude that the relevant O/H ratio in M33 is between 0.1 and 
0.4 dex below that of the Milky Way. This difference is much smaller than that between 
the SMC and the Milky Way, for which log(0/H) differs by ~ 0.8 dex. If the relatively 
small difference between the metallicity in the Milky wand M33 is responsible for the lack of 
methanol masers in the latter, it suggests that the maser luminosity must be a very senstive 
function of the galactic metallicity. 



Henkel et al. I (119871 ) have detected thermal methanol emission from two galaxies, IC342 
and NGC253, finding that the fractional abundance of methanol is ~10~ 7,5 . This is similar to 
that found in GMCs in the Milky Way, but is a factor of at least 100 lower than that required 



for hig h brightness methanol maser luminosity as discussed by ISobolev. Cragg. fc Godfrey 



( 119971 ). However, this difference is also found in comparing hot cores in the Milky Way, 
presumed to be the sites of methanol masers, with more extended molecular cloud regions. 
The methanol abundance may be greatly increased in regions near hot stars by thermal 
desorption of molecules frozen onto grain surfaces. We have no direct evidence regarding 
the methanol abundance in M33, so it is possible that the lower O/H ratio does yield an 
insufficient methanol abundance to produce highly luminous masers. 
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There are ind ividual regions within M33, presumably powered by massive young stars 
( IHinz et al. 1120041 ). which are powerful far-infrared sou rces. The infrared fl u x is a critical re- 
uirement f or maser pumping as elaborated in m odels of lCragg et al. I (119921 ) , ISobolev &: Deguchi 
19941 ). and ISobolev. Cragg. fc Godfrey I ( 19971 ). The transition to the second torsional ex- 
cited state occurs at a wavelength of ~ 30 / mi, so that dust temperatures of at l east 150 K are 
required to achieve high maser brightness (ISobolev. Cragg. fc Godfrey 1 119970 . Among our 
targets was number 8 of lEngargiola et al.l (120031 ). whi ch lies within 15"o f the optical nebula 
NGC 604. It has an IRAS luminosity of 6.8xl0 7 L^ dRice et al. 111990) and the associated 
GMC has a mass derived from CO of 4.4xl0 5 M ( lEngargiola et al.l 120031 ). By Galactic 
standards this region would seem likely to harbour a high luminosity methanol maser, but 
no such emisson was detected. 



Fix fc Mutel I (119851 ) were unsuccessful in a search using the VLA for highly luminous 
OH masers in M33. Their limit of ~ 25 mJy (5a) was sufficient to elminate the presence of 
any type I maser as luminous as the brightest type I masers in the Milky Way The present 
work thus reemphasizes the mystery of the lack of luminous masers in M33. 

We are grateful to Phil Perillat for many discussions regarding the performance of the 
Arecibo telescope, and routines for supplying routines for data reduction. We thank German 
Cortes-Medellin and Lynn Baker for their contributions to the receiver system development, 
Kurt Kabelac and David Overbaugh for machining most of the dewar components, and 
Rajagopalan Ganesan and Lisa Locke for calibrating the C-band high receiver used in this 
work. We are also grateful to numerous other members of staff at the Arecibo Observatory 
who helped us in the receiver installation and calibration, and in setting up our observations. 
This work was supported in part by the Jet Propulsion Laboratory, California Institute of 
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